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FOREWORD 



(U) This report on Project 055OT** will be published In three volumes. 
Volume I contains the theory, sample application, and the associated mathe- 
matical model. Tills effort commenced on 15 June I965 and was completed on 
15 January 1966, Volume? II and III will contain a few sample comparative 
analyses and the Energy-Maneuverability Diagrams of several aircraft fur clsan 
end air-to-air configurations. This report supersedes APGC-TP3-64-35 (refer- 
ence 8), 
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Mathematical Services Laboratory, and to Hiss Butty Jo Salter, Illustrator, 
of the Graphics Section, 
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UUClASSlTlld ABSTRACT 



TUis report shows how an aircraft's energy state and energy rate capa- 
bilities arc directly related to operational maneuverability and efficiency 
in terms of energy-maneuverability theory. It demonstrates also how energy, 
maneuverability theory may be applied to assist the tactician, commander, 
planner,' and designer in optimizing aircraft performance. Load factor versus 
velocity (G-Y) and altitude versus Mach number (H-H) diagrams are employed 
to obtain the interacting energy relationships fundamental to energy- 
maneuverability theory. The G-V diagrams provide a measure of instantaneous 
maneuverability while the H-M diagrams (the Most valuable diagrams) show 
sustained {Maneuverability as 3 function of energy rate! g, efficiency, and 
range throughout an aircraft's performance envelope, the enorg/ diagrams, 
as the working tools of energy-maneuverability theory, may be used to deter- 
mine operational maneuverability and efficiency of various armament-engine- 
airfra^ie combinations. 



This document is subject to special expert controls and each transmittal 
to foreign govenv^ents or foreign nationals may be made only with prior 
approval of APGG (PG70), Eglin APB, Florida* 
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INTRODUCTION 



(U) Aircraft maneuverability can be defined as the ability to change 
direction and/or magnitude of the velocity vector* While this definition 
describes maneuverability accurately, it provides little feel for the fighter 
pilot or engineer on how to acquire best (optimum) maneuverability. However, 
from experience, we know that the beat way to maneuver for position advantage- 
or to deny this same advantage to an opponent depends on the type of ordnance 
used end the performance of the aircraft. The type of ordnance employed 
determines the possible delivery conditions needed to effectively deliver this 
ordnance, whether it be guided missiles, guns, or bombs. Quantitatively, 
these delivery conditions can be depicted by launch or firing envelopes* Once 
the initial delivery conditions are known, the problem becomes one of maneu- 
vering into tbs effective launch envelope. Such maneuverability is dependent 
upon the ability of the pilot to control turn, altitude, airspeed, and accel- 
eration, 

/ 

(U) h The purpose of the following discussion is to show how energy- 
maneuverability is related to operational naneuverability end hpw this rela- ^ 
tionship may be exploited by the tactician, commander, planner, or designer ' 
in developing valid maneuvering and/ or delivery tactics along with better 
aerial combat weapons syataos. 
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ENsmiAflEous miziNERmum 



. (U) Turn can be described in terms of radius (r) and/or rate (tu) at 
various airspeeds (V) and radial 3 (N,) by employing the relationships 

in conjunction with aerodynamic force system equations. From such equations, 
numerous charts depicting turn radius and rate can be developed to provide 
butiie measure of maneuverability, Needless to say, the numerous charts and 
associated contours are difficult to digest. Tor simplification and clarity, 
load factor versus velocity (C-V) diagrams are employed to depict turn in a 
manner conaisteni, with a pilot's background and his cockpit instrumentation, 
(See Figure 1.) 

(IJ) The intent of this diagram (Figure 1) is to enable a pilot to deter- 
mine maximum turn in terms of g or load factor by consulting the aerodynamic 
limit at the left and the structural/stabilator limits at the top, bottom, and 
to the right. By an overlay comparison of G~V diagrams, a pilot can determine 
if he, or a possible adversary, has a turn advantage, Any turn capability or 
advantage, extracted from such a diagram, provides only a relative measure of 
instantaneous maneuverability , The diagram fails to indicate the effect of 
pulling g in terms of losing or gaining altitude and/or air3peed. As 3 result, 
no measure of sustained maneuverability can be acquired from a study of this 
diagram. To develop thia information, a look in a different direction is 
necessary. 
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SUSTAINED MNEUVERADIUTY 



CCNERAL 

4 

f U) Altitude, airspeed, and changes thereto are directly dependent upon 
the force system icting along, and normal to, the flight path. Dy mathemati- 
cally manipulating the expressions describing this force system, altitude (h) 
and airspeed (V) can be combined in the expression for specific energy (E,)^ 




as shown in Appendix II. 

(U) To i..jneuver for a desired change or a combination of changes in 
direction, altitude, and airspeed, a pilot must disturb the force system sur- 
rounding his aircraft* TEierefore, from the above expression, we deduce that 
maneuverability is not only related to directional change (turn), but is also 
related to specific energy in terms of altitude and airspeed. From this ex- 
pression, ve can also deduce that all maneuvering will be conducted between a 
maximum energy level associated with a best altitude-airspeed combination and 
a minimum anergy level associated with zero altitude and minimum airspeed, 
These maximum and minimum energy levfcLs.ifliay be represented in an altitude 
versus Mach number (H-hf) diagram (Figure 2)* The maximum energy level is 
located on Figure 2 at th^ ptfint where the specific energy (E,) contour is 
tangent to the steady-state envelope* The minimum energy level is located 
on Figure 2 at sea level where the appropriate specific energy contour inters 
cepts the steady-state envelope, (The steady-state envelope* is defined as the 
level flight operating boundary determined by angle-of-attach limits , thrust 
available, drag, and structural limits,) 

(U) In an air-to-air battle, offensive maneuvering advantage will belong 
to the pilot who can enter an engagement at a higher energy level and maintain 
more energy than his opponent while locked in a maneuver and counter-maneuver 
duel* Maneuvering advantage will also belong to the pilot who enters an air- 
to-air battle at a lower energy level, but can gain more energy than bis 
opponent during the course of the battle. From a performance standpoint, such 
an advantage is clear because the pilot with the most energy has a better 
opportunity to engage or disengage at his own choosing. On the other hand, 
oncrgy-loss maneuvers ean be employed defensively to nullify an attach or to 
gain a temporary offensive maneuvering position* Implicit in the entire 
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Figure 2 # F^4c Steady-State Envelope. 



discussion on energy state and/or energy rate advantages is the fact that a 
pilot has enough internal energy (fuel) available td exploit these advantages, 

-CH). Jp.m a lr-t o-surf ace role, a pilot is not us interested in a high 
energy state as he is in maintaining energy while maneuvering with a wide 
assortment of stores on board. If he cannot maintain maneuvering energy, his 
choice of tactics/techniques becomes limited. In addition, if this same pilot 
is tapped by enemy air, his ability to evade or nullify the attack beeches 
questionable. 

( U) Observing the correlation of energy with maneuverability, it follows 
that tactical maneuverability is related to the amount of energy possessed and 
how well that energy is managed. From a design standpoint, this means a fighter 
pilot must be given a vehicle herein such factors as energy state, energy rate, 
and He quantity of internal energy available are properly considered. For best 
maneuverability, 'the fighter pilot must know when arid how to move to a higher or 
lover energy level and how tc best conserve his internal energy when locked in 
an air-to-air or air- to-surface encounter. 

5 
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ENERGY RATE 



(U) Tor an offensive maneuvering advantage, a fighter pilot must be at a 
higher energy level or be able to gain energy more quickly than his adversary 
before the maneuver and counter-maneuver portion of the battle begins. To gain 
energy more quickly — once SCI, radar, or visual contact is made — necessitates 3 
best path for accomplishing this task. 

fU) An approximate method for finding a best flight path was discovered 
by H $ r Rutovski (sec Reference 1) t Using his method, 2s outlined in 1 Appen- 
dix II, the best (Rutovski) path for gaining maneuvering energy may be repre- 
sented on an altitude versus Mach number (I UN) diagram containing energy rate 
(specific excess power) contours within the steady-state envelope, as shown In 
Figure 3, Energy gain is maximum at the points where the specific energy (E t ) 
contours are tangent to the specific excess power (P,) contours, where 1 




T t - thrust available, D = drag, V ^ velocity, and W ^ weight, A glance 
at Figure 3 shows a best (optimum) path for gaining energy most rapidly, Not 
normally shown is the best path when the starting point is located off the 
basic Rutowskl path, A solution bo this problem becomes easy if the energy 
rate, off the Rutovski path, inaide the steady— state envelope is assumed to 
be zero. Under this assumption, the pilot moves along the specific energy 
contour consistent with his energy level until intercept is made with the 
Futowski path. As shown in Figure 3> the best path consists of two segments; 
the appropriate specific energy contour and the basic Rutowski path, Using 
this procedure, pilots can determine the best paths from any point in the enve- 
lope, However, these paths are approximate for tso reasons: (1) load factor 

ia assumed constant (1 g) in developing the ba$ic Rutowski pati and (2) energy 
rate is assumed to be zero in developing the best path from an/ point in the 
envelope. 

(U) To provide a more exact solution, A* El Bryson and H* J. Kelley (Refer- 
ence Z and Appendix II) have developed « direct method while H, P, Hermann (Ref- 
erence 3) has developed an indirect method for finding best paths. Flight paths, 
determined by these methods, show that Rutowski is very nearly correct* Rutowski’s 
method, when compared with the Bryson-Kelley and Keennann methods, reveals that 
a rule-of-thumb technique can be used by a pilot or engineer to find best energy 
paths* (See Appendix III,) The technique uses the simple rule 4E, - kAH for 
finding intercept curves and the subsonic-supersonic transition curve to the 
Rutowski path, (See Figure 4.) The value of k <= 2/3 when Mach number must bo 
decreased and k n -1 v>hen Mach number must be increased to intpreept the basic 
Rutowski path. 

6 
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(U) Even though the path developed by this procedure way be satisfactory 
to the engineer, it still is not good enough to unable the pilot to fly the 
path, because of the constantly changing altitude, Mach number, and pitch angle. 
Observation and analysis of the path Just defined, however, suggest a way to 
avoid this predicament, Generally speaking, the subsonic portion of the path 
can be represented by a constant Kadi number climb, while the supersonic portion 
may be approx imatpd by an average constant calibrated airspeed. To intercept 
rhe subsonic or supersonic se^nents, the pilot pulls up or pushes over, as indi- 
cated by the rule-of-thumb, until he intercepts the basic path. At intercept, 
the pilot should lead the Koch number or calibrated airspeed to prevent a tem- 
porary loss of energy by pulling too much g. Intercepts from the subsonic 
segment to the supersonic segment of the Kutcwski path should be accomplished 
A less than 2 g, while intercepts to the subsonic portion of the Rutowshi path 
should be accomplished at less than 3 or 4 g at lover altitudes and should 
decrease correspondingly as altitude increases, 

(U) Although the If-H diagram is useful for approximating best energy rate 
flight paths, observation reveals that it can also be used for another purpose. 
The contours contained within the steady— state envelope provide a measure of 
the ability to gain energy throughout the envelope. Since gaining energy is 
related to maintaining maneuverability, the Energy Rate diagram provides a 
measure of sustained maneuverability as a function of energy rate, fly overlay 
techniques, the L-g Energy Rate diagram can be used by the fighter pilot or 
tactician to determine if he can gain energy more quickly than some adversary. 
Actual time values, depicting how rapidly the transfer takes place,- can be 
provided by the previously-mentioned optimization programs. Such values will 
be provided In 1r Tacticai Applications, Ir Section IV of this report. When this 
infuimaUoa is cgrrelated with some analysis yet to be presented, the pilot 
or tactician can then determine the type of tactics or maneuvers to employ. 
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( U) Energy Rate diagrams of more than 1 g can be helpful in determining 
the best tactics to employ in the maneuver and counter- maneuver portions of the 
fight. As shown in Figures 5 and 6, these diagrams contain both positive and 
negative energy rate (P,) contours within the steady-state envelope, As such, 
these diagrams portray the ability to maintain energy while pulling g; hence, 
they provide a measure of sustained maneuverability as a function of g. 

( U) Once again, by simple overlay or comparison techniques, regions of 
energy advantage and disadvantage can be easily determined. If a fighter pilot 
can gain energy more quickly or lose It less rapidly than some adversary In a 
maneuvering fight, ke has offensive maneuvering advantage* On the other hand, 
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if the energy values ore reversed, the pilot; although forced on the defensive, 
may employ energy loss maneuvers to his advantage. In either case, the p-g and 
5-g Energy Rate diagrams graphically portray capabilities and limitations in 
the maneuver and counter-maneuver portions of the fitfht. In the a i r-rto-.fi ur face 
role, these diagrams may be employed to determine maneuvering capabilities and 
limitations with a wide assortment of stores on board. With this information, 
pilots and tacticians can develop pre-attack and post-attack tactics and maneu^ 
vers against a hostile surface complex. 

(U) Even though the 3-g and 5-g diagrams berve as useful tools to deter- 
mine advantages end disadvantages, they do not specify the exact tactics or 
maneuvers needed. To decide what maneuver* should be employed, a pilot must 
be veil versed in the theory, and proficient in the practice, of air-to-air and 
air-to-surface tactics (see References 4 and 5), With this background, a pilot 
can translate relative energy gain or loas relationships into valuable tactical 
maneuvers. 

( IT) Recently, the Bryson-Kelley method has been employed to develop best 
three-dimensional maneuvers (see Reference 6). This method appears promising 
in finding specific optimum maneuvers fur change of direction, rate of closure, 
and combinations thereof in minimum time or with minimum fuel until weapone 
launch. However, as presently developed, this method fail3 to consider; (1) 
the best relative regions -within the flight envelope to maneuver and counter- 
maneuver against a knowi adversary; (2) plausible counter-maneuvers by an 
adversary as he observes and/or anticipates the optimum maneuvers; (3) a 
sequence of plausible counter-maneuvers or maneuvers by an adversary 'for which 
a sequence of optimum maneuvers or counter-maneuvers will be necessary; and 
(4) the possibility that more than one optimum maneuver or counter-maneuver 
can be flown against a specific counter-maneuver or maneuver, 

(U) Because of these serious deficiencies, the Bryson-Kelley method can- 
not be used by itself to develop valid tactics for the air-to-air battle, how- 
ever, there may be a possibility of developing near optimum tactics if the 
qualitative knowledge of the tactician concerning plausible maneuvers and 
counter-Baneuvers T is used in conjunction with the quantitative methods devel- 
oped by Rutowski, aryson-XeUey, and fleemianiu The Deputy for Effectiveness 
Test, Air Proving Ground Center, and the Air Force Armament Laboratory, Research 
and Technology Division (AID), Eglin Air Force Base, Florida, are investigating 
the use of these methods for this purpose. 



EFFICIENCY 

(J) Until now, the discussion has been concerned with energy state (h,V) 
and/or energy rate (p,) in an effort to describe maneuverability and to gain ma- 
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iieuvk;eiu« iulviuitnge* Knowing l(i;it energy slate and/or etlcrgy r^te advantage Jo— 
j^rnls upon the internal energy (Fuel) available, wo kill how consider the amount 
oi paternal energy that can be converted into maneuvering energy* To acquire this 
itiiunimt ioiij a iiiuihcui.'ilieal expression must be developed which considers specific 
eitt'rgy gained versus intenial energy expended. Such an expression is 
P? 

E-MC - — w r 
*r 

where E-Hi; ^ energy-maneuverability efficiency, P* = average specific excess 
putter, * fuel weight rate flow, and v r = fuel weight* (See Appendix II for 
detailed development of this expression,) Tt*i types of E-M Efficiency diagrams, 
incorporating these efficiency contours vitliin the steady estate envelope, may 
be constructed. In the first diagram (see figure 7) constant fuel weight (5o4 
internal) is assumed. In the second diagram (see Figure 8) only the fuel re- 
maining at a given energy level is considered, after reducing the quantity of 
fuel by the minimum amount of fucl^ needed to reach that energy level. The 
efficiency contours depicted on this diagram consider fuel available minus % 
internal fuel and 30 minutes fuel for best loiter speed at ID, 000 ft. Doth 
of those E-tf Efficiency diagrams cun be used to: (l) find the most effi- 
cient (minimum fuel) paths by employing the same rule-of-thumL- techniques 
used with t^ie Energy Rato diagrams and (2) determine the amount of internal 
energy tfrat can be converted into maneuvering energy as well as the effi- 
ciency of that conversion. Since the diagrams can be employed in this fash- 
ion, they provide a measure of sustained maneuverability as a function of 
efficiency. In addition, C-H Efficiency diagrams ean bo used extensively to 
determine relative advantages and disadvantages of competing transport and 
bomber designs. For these type aircraft, load factor and energy r n -e arc 
loss important measures uf operational performance, The second E Efficiency 
diagram is more meaningful, since variable fuel weight is considered through- 
out thu Flight imvulopu. However, the constant fuel E-H Efficiency diagram 
is important in determining regions of best efficiency, independent of fuel 
histories. Thu relative merit of these two diagrams will be discussed in 
"’Tactical Applications," Section IV of this report, 

\ * 

(U) Uy employing comparative techniques, the tactician can generally see 
whether n fighter pilot or his adversary will conserve a greater percentage of 
fuel in moving from one energy level to another. Higher numerical values indi- 
cate a greater percentage of fuel remaining or a smaller percentage of fuel 
consume!). Thus, by correlating the E-H Efficiency diagrams with the Energy 
Rate diagrams, the tactician can determine to what degree a pilot or his ad- 
versary can realistically maintain or employ ony energy state and/or energy 
rate advantages, To assist in this endeavor actual fuel percentage values 
can be provided (by the optimization programs) to show how efficiently the 
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energy transfer takes piece. Such Information will be provided in ,T Tactical 
Applications/* Section TV of this report. 




RPU CE 

fU) Thus far, maneuverability has been described directly as a function 
of energy rate, g* end efficient^, However* to completely describe maneuver- 
ability, ne must consider indirect as well as direct influences* Range in- 
directly influences maneuverability as it plays a vital part in determining 
the area of maneuverability available over the earths surface* Because of 
this relationship a combat pilot must have a good but simple measure of his 
available range at any altitude^airspeeil combination within the steady-state 
envelope. To gain this information* wc must consider; (1) the fiiel consumed 
and the distance traversed in reaching any altitude-airspeed combination and 
(2) the remaining range available as a function of the fuel remaining at any 
altitude-airspeed combination, 

(U) By properly considering this information, as outlined in Appendix 
II, an tf-K diagram depicting range can be developed (see figure 9 }* From 
this diagram a pilot can determine range at any altitude-airspeed combi- 
nation including the distance traversed to reach that cc nib i nation* The range 
contours depicted on this diagram are based on ‘he same fuel reserve considera- 
tions used in the variable fuel E-H Efficiency diagram , The shaded area on the 
chart represents a transient region in the flight envelope. In this region, 
aircraft drag, l,e,, thrust required, is greater than military thrust available 
and less than minimum'^ ft erburner thrust. For this reason* steady-state # 
flight is not possible unless some device* e,g.*j3peccTbraltes ; ia employed to 
increase drag, 

(U) By using the Range diagram in conjunction with the other energy- 
maneuverability diagrams* the tactician can determine to that degree a pilot 
or his adversary can realistically gain advantage consistent with distance 
from friendly airfields or tanker support. In addition, planners and de- 
signers can evaluate the true operational performance of transport and bomber 
type aircraft by considering the Range diagram along with the E-H Efficiency 
diagrams. 
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SECTION IV 

TACTICAL APPLICATIONS 



(U) If the energy-maneuverability concept is incorporated in the present 
Wnuvlctlgc on Tighter tactics, a pilot can be provided more T meaningful infor- 
matiot] on how lie uhoulif maneuver to gain advantage. The resulting information 
wlLI then reveal to the pilot how he should best exploit the maneuvering capa- 
bilities of his aircraft. Additionally, valid comparative analyses can be 
peri'uimeil if similar diagrams are constructed fni potential enemy fighters. 

The relative regions of advantage or disadvantage are found in terms of g, 
energy vale, efficiency, and range by performance comparisons throughout the 
flight envelope. Frorn this comparison; the tactician or pilot can easily 
deLermino which of two aircraft has the advantage in terms of instantaneous 
maj i cover ability, sustained rear euvar ability, and range. Using this information, 
die tactician can determine how to best maneuyer for advantage. 
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For a sample comparison, we shall consider the F-*iC versos the Soviet 
MIG-21 and determine maneuvering advantages and disadvantages. The conditions 
for comparison will Ge typical air-to-air configurations, with 5 d£ internal 
fuel, unlessi^ecified otherwise. In the G-V diagrams { Figures 10 and 11), the 
aerodynamic g limit oT the MIG-21 lies to the left of the same limit "for the 
F-hC. At a glance. Figures 1A and 11 indicate the MIC-21 has an enorTroi^s^ffC 
stantaneous maneuverability advantage over the F-l+C, These diagrams^also in- T 
dicate that the HIG-21 has an advantage when comparing structural limits. The 
1-g Energy Rate diagrams ( Figures 12 and 13) show that^the HIG-21 has the ad- 
vantage within most of the subsonic portion of the flight envelope and through- 
out all of the supersonic portion of the flight envelope. The only region of 
advantage for the F-^C lies in the subsonic 'ml transonic areas below 13,000 
Teet, The magnitude of the maximum power energy rate advantage can be deter- 
mined by. consulting Table I, * 



TABLE. I. \ 



RUTCWSKI MINIMUM TIME PATHS (MAXIMUM POWER) 
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Figure lo. F-^C G-V Diagram at 30,000 Feet* 







Figure II. 



WIG-21 G-V Diagram at 30,000 Feet. 
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In Figures through 17, ve note that the T-kC retains most of its lov-altitude 
subsonic/transonic sustained maneuvering advantage as g is increased from 1 to 5. 
In addition, those figures reveal that the MIC-21 not only has a supersonic ad- 
vantage, but also is gaining regions of advantage subsonically since it'ean pull 
g in regions where the F-kC cannot operate, Essen tially, such a condition indi- 
cates the NIG-21 can turn more quickly than, or inside, the F-kC, This same 
conclusion was readied by studying the <LV diagram. Logically, this means the 
left-hand boundaries of the and 5-g Energy Rate diagrams provide a 

measure nf instantaneous maneuverability, while entire diagrams provide a 
measure of sustained maneuverability as a function of energy rate and g„ 

l^S) The military paver Ug Energy Rate diagrams (Figures 18 and 19) show 
that the F-kc has the ability to gain energy more rapidly than the MIG-21 
throughout most of the envelope. The magnitude of this advantage is depicted 
in Table XI. 



TABLE II, RUTOK'SKI MINIMUM TIME PATHS (MILITARY POWER) 



' 


E, = 3,000 ft to E, = 45,000 ft 


Type 

Aircraft 


Time 

(sec) 


Fuel 

Weight Used 

(lb) 


Fuel 

Used 

(» 


f-Ag 


509 


1341 


12 


" ^-tuG-ai 

flrpTT 
VJ- BmT 


403 


645 


17' . 



"(■S) By consulting the military power 3-g and 5-g Energy Kate diagrams 
(Figures 20 through 2J), we observe that the F-kC has a sustained maneuvering 
advantage at the lower altitudes and higher Mach numbers. These diagrams also 
reveal that the HIG«21 regions of advantage spread to the lower portion of the 
envelope as g is increased from 1 to 9, In addition* these diagrams (Figures 
20 through £3) reveal that the MIG-21 can maneuver in regions unavailable to 
the F-kC, As mentioned previously, such a condition indicates that the HIC-21 
can outturn the F-kC. 

From Figures £4 and £5, we observe that the F -bG has the subsonic 
advantage, while the HIG-21 has the supersonic advantage in terms of the 
maneuvering energy gained versus the percentage of internal energy expended. 

On the other hand, the variable fuel E-M Efficiency diagrams (Figures 26 
and 27) reveal that the F-kC increases its subsonic advantage and acquires 
an advantage through most of the supersonic portion of the envelope. A 
natural question arises at this point as to why the F-kC appears to have 
a greater degree of advantage in the variable fuel diagram than in the 
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Figure L6, F-4C Maximum Power 5-G Energy Rate Diagram, 




Figure 17. H1G-21 Maximum Power £-G Energy Rate Diagram* 
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Figure 22, F-^C Military Power 5-G Energy Retc Diagram, 




Figure 25* HIG-21 Military Power 5-G Energy Hate Diagram, 
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Figure 261 Maximum Power Variable Fuel E— )■( Efficiency Diagram* 
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constant fuel diagram. The greater advantage attributed to the F^C,in the 
variable fuel diagram results from the fact that the has expended a 
smaller percentage of its.^uel when it readies the regions wherp^the MI G-21 
is more efficient. Certainly such an advantage would not be realistic if the 
F-hC entered the engagement with a smaller percentage of fueJ on board and/or * 
the engagement started at a high subsonic energy level. If such a condition 
existed^ the efficiency advantage attributed to the HIC-21 in the supersonic 
region of the constant fuel diagram ^ould be realistic. This point of dif- 
ference becomes important to the F-ic pilot he encounters a H1G-21 in 
hostile territory. From the military power E->j Efficiency diagrams depicted 
\ io Figures 20 through 31, we see that the F-4C has the advantage in both cases, 
assuming, of course, that each aircraft has the same assumed or starting fuel 
percentages on board, Actual time and percentage values of fuel expended 
along minimum fuel paths between energy levels for these two aircraft are dis- 
played in Tables HI and IV. 



TABI£ III. RUIOWSKI MINIMUM TUEL PATHS (MAXIMUM POWER) 
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'TABLE IV. RLfTOWSKI MINIMUM FUEL PATHS {HIUTARY POWER) 
~T ~ l.L 3 . 000 " ft toT. = 1*5.000 ft 




Type 

Aircraft 


Time 

(sec) 


Fir el 

Weight Used 
(lb) 


Fuel 

Used 

(*) 


, F-W 


325 


1292 


u 


\ HIG-21 


1)22 

1 


6lfl 


16 



* '(■'£) The Range diagrams (Figures 32 and 33) reveal that the FJn; has a 
substantial advantage in the subsonic portion of the envelope and a lesser 
advantage in the supersonic portion of the envelope. 
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Figure 50. F-4C Military Power Variable Fuel C-H Efficiency Diagram. 




Figure 51, HIC-21 Military Power Variable Fuel E-M Efficiency 
Diagram. 
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By comparing the rule-of-thumb performance of the AIM-9B/AA-2 and 
the AIM-7E (see Appendix I), *e find the F-4C with four AIH-7E missiles has 
an enormous all-aspect, first-shot advantage aver the HIG-21 equipped with 
internal gun(s) and two AiL-2 missiles. This advantage prevails against either 
3 maneuvering or nonmaneuvering HIC-21. , However, in spite of this advantage, 
an F-bC pilot my find it difficult to employ AIH-7E du -ing a re-attack or in 
an effort tc nullify an attack, since the HIG-21 can easily outturn the F-hC 
as well as maintain more energy tfiile doing bo p By exploiting this dual ad- 
vantage, a skillful HIG-21 pilot may prevent a successful A1H-7E missile 
launch by simply maneuvering away from the front toward the rear hemisphere 
of the F-bC. For close-in maneuvering, the F-bC can mount a £0 millimeter 
centerline* gun pod in addition tc the four AIM-7 r missiles in an effort to get 
inside the missile minimum , firing range restrictions* r However, such a fix 
results in an aven greater margin of maneuvering superiority for the HIG-21 
by reducing, the already inferior instantaneous and sustained maneuverability 
of the F~4C. The magnitude of this maneuverability loss for the F-JfC can be 
determined by consulting the energy diagrams in Volume III of this report, 
to be published at a later date. 

FTom the foregoing analysis, it is clear that the HIG-21 enjoys an 
enormous instantaneous maneuvering advantage and a substantial sustained 
maneuvering advantage in terms of energy rate and g throughout the supersonic 
portion of the flight envelope - Subsonically, at both maximum end military 
power, the HIG-B1 has a sustained maneuvering advantage in the upper portions 
of the envelope that spread to the lower portions as g increases. Cn the 
other hand, the F-**C has a sustained maneuvering advantage in terms of effi- 
ciency throughout the entire subsonic portion of the envelope extending 
through most of the supersonic envelope- Only in range and first-shot capa- 
bility does the F-4C enjoy a substantial advantage over the HIG-21, 

(U) Naturally, for a complete analysis, additional in format ion must be 
developed. The tactician needs energy-maneuverability diagrams for various 
type combat configurations. In addition, he needs comparative missile fixing 
envelopes together with radar and maneuvering constraints that may be imposed 
on the pilot, or radar operator* If this information is provided, the tacti- 
cian can design tactics by uaing energy-maneuverability methods, Assuming 
that Foreign Technology Divisions can provide reasonably accurate data con- 
cerning enemy perfcflnaance, the tactician, for the first time, can develop 
effective tactics against any adversary. In addition, tactical commanders 
can use the energy-maneuverability comparative analyses to gain meaningful 
perspective for decisions concerning the employment of friendly fighters 
against a known enemy. 
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SECTION V 
REQUIREMENTS 



(U) Presently, in orddr to meet mission requirements, Air Force planners ■ 
direct that new designs meet certain specifications in terms of altitude, air- 
speed, acceleration, g, and range. Contractors, in an effort to satisfy the 
customer, produce designs to meet these specifications, However, no guarantee 
can be made that the design selected will be the best one since such specifi- 
cations are p*int data (derivatives) and provide no indication of an aircrafts 
integrated performance and design efficiency throughout the flight envelope, 

(U) However, by applying energy-maneuverability techniques, along with 
other information deemed necessary by the tactician, planners would have the 
advantage of looking at complete performance (including the previously men- 
tioned point data) before making decisions concerning aircraft requirements, 

As a result, true operational need would be considered by both planners and 
designers in determining the best overall combination of armament, engine, and 
airframe in future designs. 
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APPENDIX I 




RULE-OF-THUMB, PERFORMANCE CAPABILITY OF AIH-7E AND AIH-JJB/AA-E MISSILES 



TO The Air Force Armament laboratory (AT AIL) , Research and Technology 
Division, AFSC, Bglin Aix Force Base, Florida, A1H-9B six-degree-of- freedom, 
digital co mputer program was employed to prnihice nvpr nnp hundred lann^ 

,_ en yp jo n ag j I E.Q. 13526' scctitm 3.3{h){4) 1 

| J in^Kaytneon Company. Bedford. Massachusetts. bv means of an 

analog computer, produced launch envelopes and a parametric study of AIK-7£ 
jnis^lea^against targe -3 also pulling ftW l g to~J g [see RerFreneeTT Jl ™A 
spot cheek by AFATVs AXH-7E five-degree -of- Freedom, digital computer prograili 
showed excellent agreement with the Raytheon Company results ( An analysis of 
the AIH-7B ond AIK-9 ^/A>U2 launch envelopes revealed the rule. -of- thumb per- 
formance presented in this appendix for these missiles* Typical launch enve- 
lopes used to develop the yulfe-of-tburob performance are shown in Figures I-i 
and T s~2. : 



By carefully noting trends or patterns, the rule-of-thumb perform- 
ance of the AIK-7E and AIM-9&/AA-2 missiles can be further simplified by 
tactical organizations for operational use* 



AIH-7C MISSILES 

MINIMUM RANGE. Minimum range f R c ^ a } for nose quarter (NQ), abeam 
(AB), and tail quarter (TQ) attacks: 



Altitude 

(ft) 


NQ 

‘(ft) 


Type of rtttac 
tiE 

(ft) 


k 

rq 

(ft) 


SI 


3,500 


7, too 


5,100 


10,000 


9,250 


3,100 


5, “wo 


20,000 


10,100 


3,600 


5,300 


30,000 


1L,1<00 


9,200 


6,300 


1*0,000 


12,0CO 


9,eco 


6,000 


For t iixti s into the attack, add 1,000 : 


ft to the above values. 


from the attack, subtract 1,000 ft. 







for turns away 
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Figure I-i. AIM-TE Missile Launch Envelope 
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Figure 1-2. AIM-9B/AA-2 Hisoilc Launch Envelope 
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W HAXIHIIH RANGE. Against a Co-Speed Nonroaneuvering Target . Max iraum 
range (R^ ( ) and angle-off (<) c< ) for a nose quarter, abeam* and tail quarter 
attack ag ainst a co-speed no maneuvering target: 

Subsonic Target - Mach G.9: 



Altitude 


N<y* 

C ft/<Jeg) 


AB 

(ft) 


TQ/-T 

'(ft/defl) 


iL 


80 , 000/10 


15,000 


9,000/30 


10,000 


82,000/10 


22,000 


13,000/30 


20,000 


82,000/12 


30,000 


19,000/30 


30,000 


02,000/25 


40,000 


26,000/50 


*■0,000 


82,000/40 


52,000 


35,000/50 



For a subsonic target Mach 0.5 at sea level, NQ/4 =^60,000 f* f 10°, A3 = 2^,000 
ft, and Tq/< - 18,000 ft/30®. For each additional 10,000 ft, add 5,000 ft to 
K(J, 8,000 ft to Afl, and 5,000 ft to TQ. 

Supersonic Target; 



Altitude 

(ft) 


HQ/< 

(ft/defl) 


All 

(ft) 


Tq/4 

(ft/deg) 


SL 


82,000/10 


_ 


__ 


10,000 


82,000/10 


12,000 


7,000/30 


20,000 


82,000/15 


20,000 


14,000/30 


30,000 


82,000/30 


30,000 


20,000/30 


40,000 


02,000/45 


40,000 


26,000/30 



N), /■gainst Honiraneuveririg Targets with Attacker Velocity Greater or 
Less than Target Velocity , Maximum ranges for nose quarter end tail quarter 
attacks against nonmaneuvering targets with attacker velocity greater or less 
than target velgclty (delta Mach) follow: 

Nose Quarter Attacks; 

1. Add 3,000 feet ta R„ a far each 0,1 delta Mach below 10,000 
feet when target velocity is greater than attacker velocity, 

2, Add 1,500 feet to R ,* t for each 0,1 delta Mach above 10,000 
feet when target velocity is greater than attacker velocity, 

3* Add 1,000 feet to R 1#1 for each 0,1 delta Mach when attacker 
velocity is greater than target velocity, 
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Tall Quarter Attacks.' 

1, Add 2,000 feet to R iax for each 0 t l delta rcach rate of 

closure* 

2. Subtract 5,000 reet frora fi atx for each 0.1 delta mach 

separation. 

l^L Against a Maneuvering Target , 

Nose Quarter Attacks r 

1, Reduce R aax 30,000 feet when target maneuvers away from 
the attack at 2 g. 

2. Reduce R AaX 5,000 feet/g for target maneuvers away from 
the attack with g greater than 2. 

’('CJ Maneuvers Away iron: Attack, Tail Quarter Attacks* 

Altitude Target G R( iax M ) 

Below 20,000 ft 3 2/3 R ix , 

Below 20,000 ft 5 1/2 R,.* 

Above 20 ,000 ft 3 1/2 R, BX 

Above 20,000 ft 5 1/3 JU** 

The above values do not include background clutter associated with a target 
at low altitude. 
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APPENDIX II 



HATHD1ATIGAL DERIVATIONS AND MODELS FOR DEVELOPING 
ENQIGY-MAN EUV GRAB IL1TY THEORY AND ASSOCIATED FLIGHT PATHS 

(Appendix II is unclassified in its entirety) 

In this appendix a discussion o£ the mathanatical methods employed to 
develop the Energy-Maneuverability Theory and associated flight paths in the 
altxtud&-Hach number plane will be presented. For conveniencej tba derivations 
will Le described in terms of the following computer programs viiich have beer* 
formulated to handle the computational aspects of the theory: 

Part I - B;,sic Energy-Haneuverability Computer Model 

Part II - The Bryson- Kelley Steepest Ascent Optimization Program 

Part III - Dynamic Profile Generator Program- 



\ 
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E^lSlC EKERG V-+1AW E UVrRABI LIT Y COMPUTER MODEL 



DERIVATIONS 



INSTANTANEOUS MANEUVERABILITY* For any given aircraft, iwinum load 
factor (normal acceleration) may be computed as a function of altitude End 
airspeed; 



where 



n L 




i 



= maximum normal acceleration (dimensionless) 




q = | p V s , dynamic pressure (lb/ft 2 ) 



p r atmospheric density (slugs/ ft-) 
V * true airspeed ( ft/sec) 



S - reference wing area (ft 2 ) 






maximum coefficient of lift (dimensionless) 



w => aircraft weight (lb) 

Since calibrated airspeed (CAS) is more meaningful to the pilot than true air- 
speed, the G-V diagrams (see Figure 1, page 5 ) depict maximum normal accel- 
eration versus CAS* 



SUSTAINED KANEUVERABILITy, Energy Rate , The energy (E) possessed by an 
aircraft is the sum of its potential energy (E p ) and its kinetic energy (E^). 
Ha thematically, 



E = E p + q. 

■ wh + 3 mV 2 




kO 




where 



h = altitude ( ft) 
m = aircraft mass (slugs) 

g = 32,l7*» ft/sec®, the gravitational acceleration 

The expression, £ = w Qi + 

of an aircraft at any altitude-airspeed combination, However, since the main 
interest lies in comparing aircraft with different veights at the same altitude- 
airspeed combinations, it is more meaningful to males the above expression inde- 
pendent of aircraft weight* Dividing both sides of the above expression by w 
yields 




The term E/w can be regarded as specific energy (E f ), with the result that the 
energy state of an aircraft can -now be expressed as a function of altitude and 
airspeed: 




The problem of managing energy involves controlling the rate of transfer 
between energy levels. Differentiating the above expression results in 




— Jt gives us a measure of the anergy state 



where the dot(-) indicates the derivative with respect to time, *j|"\ To 
provide more insight into energy rate, £„ , we may employ Figure 11-1 and write 
a force balance equation along the flight path. 



or 

or 



^ T l _ D - w sin y , 
v * 

T, - D - v sin y *- -r V, 
S 



T*-D . V 

— - iin v+ y 



Multiplying both sides of this expression by V yields 




W 

Vsin v+-. 
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rifjure II-l, Aircraft Force-Balance Diagram* 
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Since h = V sin y, ve may write 




a 

The right aide of the above expression is equal to E ># Recalling that 
work is accomplished in transferring frcm one energy level to another* and 
that power, by definition, is the time rate of doing work, the left side of 
the above equation may be equated to specific excess power, P, : 




In an attempt to counter an immediate threat, the energy-oriented fighter 
will strive to increase his maneuvering energy as quickly aa possible* 
This amounts to maximizing the rate of transfer between energy levels, which 
is equivalent to maximizing the integral 




According to Rutovski (reference l) f this is accomplished when 




In the altitude-Mach number plane, these relationships are satisfied ttt those 
points where the E f contours are tangent to the L*g P* contours, Connecting 
these points results in an approximate minimum time path* 



Energy-Maneuverability Efficiency (IUHE) . If the above-mentioned threat 
ia not as imminent! the pilot will attempt to increase tils maneuvering energy 
while conserving internal energy (fuel) for future maneuverability. This is 
achieved by maximizing the integral 



E '-=f.r‘ 



dw. 







Since dE. 



spd 



. - P.4t„ 

* * 

dw a dt (w f :■ fuel flow - Ib/sec) 



ve see that 

dE, P, 
rfv “ ’ Jl ' 




Again, by employing Rutowski^e technique, we obtain 




These relationships are satisfied at those points in the altitude-Nath number 
pljne where the E, contours are tangent to the 1-g P t /w f contours. Connecting 
these points results in an approximate minimum fuel path. 

The P B />^ contour e suggest a measure of efficiency In view of the fact 
that they depict the amount of specific energy gained per pound of fuel ©u. 
pended. In order to acquire a more meaningful measure of efficiency, these 
contours can be modified to portray the amount of maneuvering energy gained 
for the internal energy (fuel) expended. This ie dtffie by multiplying the 
P*/4 contours by the weight of fuel available, to obtain the resulting 
expression for Energy-Maneuver ability Efficiency; 

t-ME •= ti v, K 

P 

t 

where PJ - * the overage P, over the fuel weight interval -fc*v t *Hr 

(ft/.ec),and h» = w fa - fc - i* ,, 

* where w f B * initial fuel vwight ( lb) 

fe c fuel consumed in flying from some reference energy level to 
any given altitude -Mach number point (lb) 

^ r '=» fuel reserve (lb) 

, . 41* 
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RANGE, For any altitude-airspeed combination, available range for cruise 
condition may be expressed as 

w M 

« = 75- V*x , 

w a 

wEiert? = the average cruise fuel flow, fr e , over the fuel weight interval 
v t0 - fc s wf s ^ r (lb/sec), 

and x =r the horizontal distance traversed in flying from some reference 
energy level to any given altitude^airspeed combination- 



COWITTATIONS 

ENERGY RATE DIAGRAMS. For any given aircraft, an Energy Rate diagram may 
be constructed by dividing the alt i tud e-Ha ch number plane into a rectangular 
grid, computing energy rate fl' ( ) values at all of the points of intersection of 
the grid lines, and then connecting points of equal R,, The contour defined 
by R, - G represents the steady-state boundary of the aircraft. An aircraft 
cannot operate outside this contour without losing energy, either in the form 
of altitude, airspeed, or some combination of both» The steady-state boundary 
is further restricted on the left by the buffet boundary (obtained by connect- 
ing points where C L = and on the right by placard limits (a combina- 

tion of pressure {structural) limits am! engine temperature limits). 

Considerable insight into the effects of "pulling g within the aircrafts 
flight envelope can be gained by constructing Energy Rate diagrams of more 
than I g, These diagrams contain both positive and negative p, cor tours within 
the 1-g steady-state envelope. As such, they provide a measure of sustained 
maneuverability as a result of pulling g within the envelope* 



L*H EFFICIENCY DIAGRAMS, Computational aspects of this diagram proceed 
in the sane manner as for the Energy Rate diagrams, except that now we compute 
and connect points of equal E-ME, Two different types of C-ME diagrams are 
constructed. The first type is referred to as the path independent (constant 
fuel) E-NE diagram. Computations for all points in the envelope are based 
on 5Dtf fuel weight. Since fuel weight is held constant, the expression 



^5 



E-HE = — vi . 
t'f 




reduces to 




The diagram Is called peth-independent since the amount of fuel at the altitude- 
Mach number points where computations are made is independent of the paths re- 
quired to reach these points. 

In the second type of E-ME diagram, called the path- dependent (variable 
fuel) E-ME diagram, the amount of fuel required to reach any given altitude- 
tfach number point is subtracted from the total fuel weight before E-ME computa- 
tions are made. The assumption ia that the pilot has flown a minimum fuel path 

from some reference energy level (we use E, * 3000 feet) to the altitude- 

REF 

Had; munler point mxJer consideration. A more detailed diecusslon of this os- 
sumption will be given later in this appendix and in Appendix III. Additionally, 
the amount of fuel upon which the path-dependent E-ME computations are based is 
reduced by a suitable reserve (normally % of full Internal plua £0 minutes 
loiter at 10,000 ft), 

RANGE DIAGRAMS, Again, the computational aspects of this diagram are 
essentially the same as for the Energy Rate and E_H Efficiency diagrams. To 
compute range, the program requires, as an additional input, a partial power 
setting table, i.e,, a table of cruise fuel flow as a function of altitude, 

Mach number, and drag (thrust required). The subsonic and supersonic* por- 
tions of the envelope ore computed using partial military and partial after- 
burner pwer settings, respectively, A transient region is observed between 
the subsonic and supersonic portions of the envelope. In this region, level 
unaccelerated flight is net possible as the thrust required is greater than 
militaiy thrust available, yet less than minimum afterburner thrust available. 

For range, only a path -dependent (variable fuel) Range diagram is con- 
structed, For this diagram, the amount of fuel available at any given altituda- 
Haeh number point is reduced by the amount of fuel required to fly a minimum 
ftiel path from some reference energy level (again, we use E_ = 3000 feet) 

REF 

to the point under consideration, and by a suitable fuel Reserve (c.g., 5^£ of 
full internal plus £0 minutes loiter at 10,000 feet). The "horizontal distance 
traversed in flying the above-mentioned minimum fuel path, x, Is considered 
part of the available renge, A discussion of the method used to compute fuel 
consumed ant) horizontal distance traversed is given later in this appendix find 
in Appendix III, 
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EXTENSIONS OF THE RUTOWSKI TECHKIQUE 



Earlier In this appendix, we saw that Rutovskl calculated the location of 
the approximate minimum time and minimum fuel paths In the altitude-Hach num- 
ber plane* Rutowskl did not, however, give any measure of the time required, 
fuel consumed, or horizontal distance traversed in flying along these approxi- 
mate paths. His two basic assumptions were that (1) the path be computed for 
1 g and (2) the weight be held constant, 

Rutowskl^s method has been extended to allow an approximation of the time 
required, fuel consumed, and horizontal distance traversed In flying along the 
minimum time or minimum fuel path. A byproduct of this extension has been the 
removal of his constant weight assumption, 

Computations proceed as follows. Program inputs Include initial, incre- 
mental, ’and final values for E a and h; 

B.. UE.) V 

1 L 

(Ah) h L . 



The specific energy equation, 

. V 2 

E f n h + —• * 

2g 

is rearranged to the form 

v= [2g<c t - h>]*. 

Then for E, = E, , E fl + AE,, E, + 2AE,, , , , , E, , the following array 
is constructed; 

hi V x ^ T. x Wi 

h a V a M a T la w, 

• » • * • ( t f a « 

till i < i i i * 



C Ei \ 01 {P 'A*lx 

C I« \ Da P * 3 < P ^>» 
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Actually the array does not run all the way from h x to h * The lower 

L 

altitudes result in Mach numbers higher than the aircrafts capability, hhen 
this occurs, altitude is incremented instead of computing T 4 , w f , , P ¥ /fy, 

TSie higher altitudes result in C's greater than C f f s. This fact eliminates 

many of the lower lines of the array, Additionally, when h > E ( , the quantity 

[£g (b\ - h)]a becomes negative, eliminating linea of the array. Finally, 
other numerical techniques, beyond the scope of this appendix, Ore employed 
to reduce the size of the above arrays, thereby decreasing the computer time 
required to construct a path. 

If a minimum time path is being computed, the program selects the altitude- 
Mach number point for which \\ is maximum in the array and this point becomes a 
point on the minimum tine path. Once the line containing the maximum l J 4 Is 
selected, it is used, along with a similar line on the previous array, to ap- 
proximate a time increment, it, a fuel increment, aw, and a 
increment, Ax, in th* following manner: 

AE / — 

it = — (the bar denotes the average, e,g t , P, 

P. ^ 

* tL ' 
bM = ■ 

(p./m 

ah ^ i 

» » [* - («)*] a * 

The method outlined above results in altitudc-Mach number points through 
which a minimum time path may be drawn* The it’s, Aw’s, and Ax’s are summed 
'over the path to give an approximation of the time required, fuel consume*!, 
and horizontal distance traversed. 

Tacit time F, is incremented, the weight used to compute the quantities in 
the array is first decremented by the quantity aw, computet! in the previous 
array. This results in a variable weight path, 

Computations for a minimum fuel path proceed in exactly the same manner, 
except that, for the minimum fuel psth, the lino is chosen in the arrays where 
p.M instead nf P f , is maximum. 

To compute a path-dependent E-MC or Range diagram, a Rutowski minimum fuel 
path must be computed first and the following table built; 



horizontal distance 



p *. * Va 
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i 

fc, . £ AHj , 

** = I “j • 

J=1 

Then, when E-ME or range is computed for a given altltude-bach number point, 

(h, N), that (h, H) determines an E„ which, in turn, determines an fc and an x 
if we assume that the fuel consumed and horizontal distance traversed In flying 
to any point on a constant specific energy line 13 the seme. This is a rather 
bold assumption, however, and cannot be accepted without further discussion. 
Appendix III provides a detailed treatment of this assumption 



PART II. Till! BRYSON-KELLEY STEEPEST ASCENT OPTIMISATION PROGRAM 



The second computer program employed in these analyses is the Brysan-Kelley 
Steepest Accent Optimization 1'rogram which provides dynamic flight profiles, in 
miniiiiunj time or with minimum expenditure of fuel, for transfer between two 
enerjjr' levels (hi, Hi) and (h 3 ,)fc). 

By using tile methods of E ,S * Rutowski, rt>n stalled Rutowski path, explained 
previously under the Energy-Maneuverability program, is obtained and is depicted 
in Figure page 7 of this report f However, these methods provide nothing more 
than very good approximations to the solution of the minimum time or minimum fuel 
problem, They provide no insight into such parameters as load fart or or pitch 
angle along the path, In addition, the? methods are predicated on 1-g level flight 
parameters and do not consider the forces acting perpendicular to the flight path* 
In essence, the Rutowski method is a static method, in itself, but a very valuable 
tool leading to the more accurate dynamic optimum paths, 

Even the more sophisticated fynainio Profile Generator Program, discussed in 
Tort III of this Appendix, provides only approximations to the desired solution 
to the minimum time or fuel problem. Admittedly, the results of using the 
Rutowski paths in conjunction with this program arc much more realistic, as now 
both load factor and pitch angle are considered throughout the path. However, 
the techniques embodied in this program are still limited by tbc alt i+ud e-Ha ch 
number combinations input into the program as points describing the approximate 
path, and yield nothing but a better approximation to the desired optimum path. 

The program is invaluable, though, as a generator of load factor an a function 
of time for input into the Bryson-Kelley Steepest Ascent frograin as the nominal 
path. 



This st eepost-as cent method of optimization is an iterative scheme which 
begins with any non- optical path and proceeds to derive a slight improvement 
each iteration frem this nominal path. This slightly improved path at each 
iteration is used as the new nauinal path, and the process is repeated until, wc 
arc sufficiently close to the optimum for our purpose. In this process, each 
new path is found by taking the trajectory which yields the largest gain in 
performance for a given size of perturbation in the control variables. 

Tlie value of a good first guess at the nominal path is immediately evident, 
If this path is close to tbc optimum, the number of Iterations necessary to 
arrive at this profile will be small indeed when compared to those necessary if 
the nominal path is far from the optimum, The ability to input good nominal 
paths results in tremendous savings of valuable computer time. 
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The analysis of single stage trajectories by the steepest ascent method 
has been thoroughly treated in the literature. One of the clearest treatments 
available is that of Bryson and Denham in Reference £* Tor convenience, a 
brief description of the general problem, as formulated by them, is repeated 
here. Some of the detailed derivations which are omitted here are presented 
in Reference 2, 

After presentation of the general problem, the specific applications 
made in formulating the program at the Air Proving Ground Center are given in 
detail. 



GENiXU PROBLEM 



Determine in the time interval ti £ t ^ t? so as to maximise 

(1) i = t # ] 

subject to the constraints 

(2) f « tfSfr), t s j * 0 

CO j£ = ftset), 3(t), tj 

£^) the given initial conditions x(ti) 

(5) t* determined by U = uLxCt^Jjta j| * 0. 



The " over the symbols above indicates a matrix quantity,, and a more 
detailed description of the above quantities fellows; 



(6) B(t) = 



(7) x(t) = 



°i(t) 

<*00 



«w(t) 

xi(t) 



x*^) 



, an in y 1 matrix of control variable programs, 
which we are free to choose. 



, <m n & l matrix of state variable programs, 
resulting from the choice of S(t) and x{ti), 
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(fi) 



■tv 

U 




a p x 1 matrix of terminal constraint functions, each 
a — function of x(t 3 } end t a , 



C9) $ - the pay-off function, a kncvn function of x(ta) and t 3 , 



■Can) * = 




an ti x 1 matrix of known functions of x(t), 5(t), and t, 
and 



fa 



(11) ft = 0 is the stopping condition that 
and is a knevn function of E(t 5 ) 



determines final time t 3/ 
and t 3 . 



The method proceeds as follows; 

1« Choose a reasonable nominal control variable program, and use 

it with the initial conditions (4) and the differential equations (3) to calcu- 
late, by numerical methods, the state variable programs x (t) until ft = 0. In 
general, this naninal pHth will not satisfy the terminal conditions T = 0, or 
yield the maximum possible value of I, 

2, Consider snail perturbations 6o( t) about the nominal control variable 
program, a (t), where 

(12) * 5(t) - 



As a result of these perturbations, the state variable program undergo 
perturbations 5x(t)j where 

m = set) - x*'(t). 

If the relations (12) end (J3) are substituted IiTto t!ie differential 
equations, given by (}), the linear differential equations 

(l4) Jj- (Ax) ^ f(t)6x + C(t)6a 

are obtained, accurate to fir?t order in the perturbations, where 





The symbol ( )* Indicates that the enclosed partial 
ated along the noninel path. 



Using the theory of odjoint differential equations, 
a ion a nay be written, 

( 17 ) d* « f *X t '(t)G(t)4a(t)dt + ^(ti)6x(tj) 

(10) d? * f a 'y(t)0(t)65(t)ut ♦ A*(ti)«(tO 

Jl i 

( 19 ) dn = J^ 3 ^(t)c(t)6S(t)dt + ^(tiJdSCti) 




‘derivatives ere evalu 
the following e^epres- 
+ Idtj 

+■ t<ft # 
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where the symbol * indicates the transpose of the matrix and elements of the 
three A matrices, appearing above, are obtained through the numerical integra- 
tion of the difftrential equations adjoint to equations ( 3 ): 

(20) 

dt 

with the boundary conditions 
where 

(a21 2! - I®!. 111 

( J 3* ' la*,.' 3x a ' S*]' 

flit ®i* 

jfixi Bxj 3x, 

|*+a }lz .t.j Bfa 

3*1 3xa dx„ 



2*2 222 ... 22» 

[Bid 3 x s 3x*| 

an fin so aril 

« " [a*/ a*] 

and 

" Hi-i’l: '-(S’S'L- 

••e-s')u- 

Hole that the Me are influence functions in that they tell hew much a 
certain terminal condition ie changed by a small change in aome initial state 
variable. Note also that the adjoint equations (SO) must be integrated back- 
wards since the boundary conditions (21) are given at the terminal point. 
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For eteepest ascent the 6o(t) program that maximizes the d* In expression 
(17) must be found, given values of dy and do ■ 0 in expressions (18) and (19), 
respectively. This maximization must also be subject to a given value of the 
integral 



(2*0 (dP)= = jS?'(t)i?(t)63(t)dt. 



The value of (dP) 2 is chosen such that the perturbations will be small 
enough to insure that the neglect of second and higher order perturbations 
leading to equation (lh) is reasonable* In addition, values of df are selected 
to bring the next solution closer to the desired terminal constraints, f rO. 

The men matrix 0(t) is symmetric and contains weighting functions as 
elements. They may be chosen arbitrarily to improve convergence* In the usual 
case (the APGC program falls into this category), k(t) is taken equal to the 
identity matrix and (dp) 2 becomes the integral of the square of the control 
variable perturbations, 65[t)* t Observation reyeala that all control variables 
should have the same dimenaioijh for equation (24) to have any meaning. To meet 
this requirement the control variables are normally required to be nondimenaional, 

A rather involved aeries of mathematical manipulations (presented in an 
orderly and clear fashion in Reference 2, but emitted here for the sake of 
brevity) leads to the following proper choice of Gp(t): 



, , \ iC^) 2 - dP'fcJdp 

M 4W :^|b M W' 



(2£) dS = df - X'ct^s^to, 



(S7) f— ~ 

H(ta; 

¥ 

(28) Ayn - Ay - Aft, 

(29) iyy - g^ndt, 
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(J°) J 4 */n S W -1 0 

(31) 

and the ( y 1 indicates the inverse matrix arid the + or - sign before the 
radical in (2J) is chosen if $ is to be increased or decreased, respectively. 

If the selected d7 is such that dp is too large> th(?n the numerator in the 
radical in (25) might become negative and a limit to the fU2e of dfJ for a given 
dP roust bejuncosed, Since dP is chosen to insure valid linearization the 
selected df must also be limited. 

The predicted change in ♦ far the change in o(t) given by (25) 

(32) « « i^(dP) 2 - <ip'lt|de)[l„ - 

+ ^jo(tl)ox(tl)* 

If df = 0 (the terminal constraints having been satisfied) and 6x(ti) = Q } 
then d? a C and equation (}2) becomes 

° 3) ® = *y i **~ 

which is a gradient in function space, since dP is the length of the stop in the 
control variable program. As the optimum program is approached and the terminal 
constraints are metj (df = 0) 3 thiB gradient must tend to zero, and expression 
(52) becomes 

(3*0 « « 

3# A nev control variable program is now obtained as 

^(35) 3ft) - 5*ft) + 

This new 5(t} is row used in the original nonlinear differential equation 
given by (3)/ and the process is repeated until the terminal constraints (2) are 
met and the gradient (33) becomes nearly zero, 
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